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We report the first search for supersymmetric particles via s-channel production and decay of 
smuons or muon sneutrinos at hadronic colliders. The data for the two-muon and two-jets final 
states were collected by the D0 experiment, and correspond to an integrated luminosity of 94±5 
pb _1 . Assuming that i?-parity is violated via the single coupling A 2 n, the number of candidate 
events is in agreement with expectation from the standard model. Exclusion contours are given in 
the (mo,mi/ 2 ) and (m^,m^) planes for A2n=0.09, 0.08 and 0.07. 
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Events with at least two muons and two hadronic jets 
in pp collisions provide a good sample in which to search 
for new physics because the contribution from standard- 
model processes to such states is rather small. Any excess 
in such topologies can be attributed to a signal from /im- 
parity violating supersymmetry (SUSY), where i?-parity 
is not conserved either in the production or in the decay 
of sparticles. 

7?-parity of any particle B is defined as R p = 
(-l) 3S+i+25 , where B, L and S are the baryon, lep- 
ton and spin quantum numbers. R p equals +1 for SM 
particles and —1 for supersymmetric partners. The con- 
servation of i?-parity is often assumed in experimental 
searches, because, without that, simultaneous lepton and 
baryon number violation would lead to rapid proton de- 
cay. However, this argument can be circumvented if lep- 
ton and baryon number conservation are treated inde- 
pendently. 

In supersymmetry, i?-parity violation (R\p) can occur 
through terms in the superpotential, that are trilinear in 
quark and lepton superfields Jl| : 

A, ,;. /., /. , I:',, + \', ,,,/.,() , I )\. + \'{ jk U?Dpl, (1) 

where i,j, and k are family indices; L and Q are the 
SU(2)-doublet lepton and quark superfields; E, U, and 
D are the singlet-lepton, up-quark, and down-quark su- 
perfields, respectively. 

Such^lp couplings offer the possibility of producing sin- 
gle supersymmetric particles , which is not the case for 
i? p -conserving supersymmetric models, in which particles 
and sparticles are always produced in pairs. Although the 
IjLp coupling constants are severely constrained by low- 
energy experimental bounds IjJQ], s-channel production 
and decay of sparticles can have a substantial cross sec- 
tion at lepton and hadron colliders ||[|. 

At pp colliders, either a sneutrino {v) or a charged 
slepton (I) can be produced in the s-channel via A- jfc cou- 
pling. In most SUSY models, the slepton has two possible 
.Rp-conserving gauge decays: either into a chargino x ± or 
a neutralino x°. These are favored over$ p decay because 
of the small value of the coupling for the latter || . Conse- 
quently, for a single dominant A^ fc coupling, production 
of a slepton (smuon or muon sneutrino) provides either a 
chargino or a neutralino, together with either a charged 
lepton or a neutrino, in the final state. 

In this Letter, we consider the resonant production of a 
muon sneutrino or a smuon via A^n coupling which leads 
to a final state containing at least two muons and two 
jets. From low-energy measurements the X' 211 coupling 
is constrained to be less than 0.06/100(GeV/c 2 ) x m,£ 
0, where m d ~ is the mass of the dfj-squark. The light- 
est supersymmetric particle (LSP) is assumed to be the 
lightest neutralino. We also assume that all sparticles 
cascade-decay into the neutralino, which decays through 
the dominant IjL p A 211 coupling. Hence, ultimately, all 



SUSY particles decay into the lightest neutralino, which 
decays into two jets and a muon. The decay of the muon 
sneutrino into a muon and a chargino, and of the smuon 
into a muon and a neutralino, therefore lead to at least 
two muons and two jets in the final state. The decay of 
the smuon into a neutrino and a chargino can also lead 
to the same topology, but only when the chargino decays 
into muon+A, and for this reason the contribution of 
that channel is small (less than 5% of the signal) and ne- 
glected in our analysis. The decay of the sneutrino into 
a neutrino and a neutralino yields only one muon in the 
final state. 

Our framework is the so-called minimal supergravity 
model (mSUGRA), which assumes the existence of a 
grand unified gauge theory and family-universal bound- 
ary conditions on the supersymmetry breaking parame- 
ters. We choose the following five parameters that com- 
pletely define the model: mo, the universal scalar mass 
at the unification scale Mx] mi/2, the universal gaugino 
mass at Mx\ A = A t = Af, = A T , the trilinear Yukawa 
coupling at Mx, sign(p), the sign of the Higgsino mix- 
ing parameter; tan (3 = < H u > / < Hd > where < H u > 
and < Hd > denote the vacuum expectation values of the 
two Higgs fields. The dependence of the cross section on 
different SUSY parameters can be found in Ref. ||. 

The data for this analysis were collected during the 
1994-1995 Fermilab Tevatron running, at a center-of- 
mass energy of 1.8 TeV, and correspond to an integrated 
luminosity of 94±5 pb _1 . The D0 detector is described 
elsewhere Q. Here, we outline the performance of the 
components relevant to this analysis. Jets are identi- 
fied using the energy deposited in the calorimeter, and 
reconstructed with a cone algorithm in pseudorapidity 
(rj) and azimuthal angle (</>) using a radius of 0.5. The 
calorimeter covers the region of |r/| < 4.0, and provides 
a resolution for electrons and single hadrons (a{E)/E) 
of 15%/a/E and 50%/v^, respectively. Muons are de- 
tected using both tracking chambers (three layers of pro- 
portional drift tubes (|?7| < 1.7), one in front of, and two 
behind magnetized iron toroids) and through ionization 
deposited in the calorimeter. The muon momentum res- 
olution is er(l/p) = 0.18(p - 2)/p 2 + 0.003 (with p in 
GeV/c). 

Events are required to satisfy a y, + jet or + jet 
trigger. The trigger efficiency is 71% and 50% for central 
and forward muons, respectively. Muons are required 
to have a transverse momentum greater than 8 GeV/c, 
and jets are required to have transverse energy exceed- 
ing 15 GeV. We apply additional criteria to select two 
isolated muons and to eliminate cosmic-ray muons. If 
there are more than two isolated muons (which happens 
only rarely), only the two leading muons are used in the 
ensuing analysis. 

The signal topologies were generated with the SUSY- 
GEN Monte Carlo program || using the cross sections 
computed in Refs. |q|| for a wide range of (too, W1/2) 
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masses. For illustration purposes, we choose a refer- 
ence point in the mSUGRA parameter space: mo=200 
GeV/c 2 , TOi /2 =243 GeV/c 2 , tan/3 = 2, A = 0, and a 
negative sign for /i. These parameters predict the fol- 
lowing sparticle masses: ma =263 GeV/c 2 , =269 
GeV/c 2 , m-± =207 GeV/c 2 , and m^o =102 GeV/c 2 . 
For A' = 0.09, the production cross sections are 1.22 pb 
and 3.34 pb for v and fl production, respectively. 

The dominant backgrounds are from tt, WW + jets and 
Z +2 jets events. The background was generated using 
PYTHIA (To), with a cross section of 5.9 ± 1.7 pb 0, 
the Z+2 jets background with VECBOS, ||| interfaced 
with the IS A JET fragmentation code p3| , and a cross 
section of 9.7 ± 0.9 pb. The VTIU±jets background was 
generated using PYTHIA Jl(| ; it provides a much smaller 
background than the tt and Z+2 jets channels. The sim- 
ulation of the detector was performed using both a full 
and a parameterized simulation. 

We use a neural network to discriminate signal from 
background in our analysis (lj] and we cross-check this 
with a more standard sequential analysis at several points 
of the SUSY parameter space. The following quantities 
are used as inputs to the neural network: the scalar sum 
of the transverse energies of the two leading jets, the 
scalar sum of the transverse momenta of the two lead- 
ing muons, the distance in (77, <j>) space between the two 
muons, the dimuon mass, the (f>) distance between 
the most energetic muon and its nearest jet, the apla- 
narity and the sphericity of the two leading muons and 
two leading jets in the laboratory frame fi5fl . 

The output of the neural network is obtained sepa- 
rately for the sneutrino and the smuon channels. The 
signal-over-background ratio for the neural network is op- 
timal for an output cutoff of 0.0 for the v and —0.10 for 
the /jl analysis. 

For the reference point, 6.42 ± 0.06 v and jl events are 
expected. The estimated background of 1.01±0.02 events 
is consistent with the two events observed in data. The 
details of the background estimate are given in Table I, 
with the quoted uncertainties being only statistical. 

The systematic errors are shown in Table || The un- 
certainties due to jet energy scale and the measurement 
of the muon pt are deduced by varying the jet Et and 
muon pt by one standard deviation. We use a fast version 
of the detector simulation for most of the SUSY points, 
and the systematic error associated with this procedure 
is also given in Table |l| The last three lines give the 
final results for the number of events, the overall sta- 
tistical error, and the overall systematic error. Using a 
Bayesian method to calculate the level of exclusion [^6| , 
our specific reference point is rejected at the 97.7% C.L. 
for X' 211 = 0.09. 





6.93 ± 0.1)5 




z.4y ± U.U4 


Total expected signal 


bA2 ± U.Uo 


tt 


0.27 ±0.01 


Z±2jets 


0.73 ±0.02 


WW + jets 


0.01 ±0.00 


Total background 


1.01 ±0.02 


Data 


2 


CL 


97.7% 



TABLE I. Number of events (expected) for the reference 
point for signal at A' = 0.09, for the background, and the 
number observed in the data after making all selections. 



Source 


Signal 


tt 


Z+2 jets 


Jet energy scale 


2% 


4% 


5% 


High efficiency 


1% 


7% 


4% 


Cross section 


10% 


30% 


10% 


Trigger simul. 


5% 


5% 


5% 


Luminosity 


5% 


5% 


5% 


fast / full simul. 


1% 


1% 


1% 


Total events 


6.42 


0.27 


0.73 


Overall statistics 


± 0.06 


± 0.01 


± 0.02 


Overall systematics 


± 0.80 


± 0.09 


± 0.10 



TABLE II. Systematic uncertainties on signal and back- 



ground, and the number of expected events, with their statis- 
tical and systematic errors. 
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To set exclusion contours, we scan the (m , mi/2) 
plane for three values of the coupling constant X' 211 = 
0.09, 0.08, 0.07, two values of tan/3 = 2, 5, all for 
sign(/i) = -1. For A' > 0.09, the coupling lOOA'/m 
is almost completely excluded by earlier experiments 
in our domain of sensitivity in m d - . The resulting ex- 
clusion contours at the 95% C.L. are shown in Figs. |l| 
and |2| in the (mo, TO1/2) plane. The most interesting fea- 
ture, is the exclusion of mx/2 values up to 260 GeV/c 2 
for tan/3=2 and A' 211 = 0.09, and the exclusion of v and 
jl with masses up to 280 GeV/c 2 . 

For low values of mo and m^, the smuon mass is 
close to the chargino or neutralino mass, the px spectrum 
of the muons is soft, and the search is inefficient. For 
(i >0 and higher values of tan/3, the sensitivity of our 
reach is expected to decrease due to the fact that the 
photino component of the LSP becomes small, resulting 
in the decrease of the branching fraction of the LSP into 
muons. In addition, charginos and neutralinos become 
light, resulting in events with softer muons and jets that 
fail the kinematic requirements. 

To conclude, a search for single smuon and single muon 
sneutrino production in the mSUGRA model with R- 
parity violation, has been performed for the first time at 
the Tevatron. We exclude mi/ 2 values up to 260 GeV 
(the excluded value of mi/ 2 depends on the value of mo) 
and sneutrino and smuon masses up to 280 GeV. The 
excluded domain in the (mo, my 2 ) plane extends the 
region excluded using the dielectron channel jl7|] . 
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FIG. 1. Exclusion contours at the 95% C.L. in the (mo, 
m i/2) plane. The top figure shows the exclusion contours 
for tan/3 = 2, \' 211 =0.09, 0.08 and 0.07. The bottom fig- 
ure shows the exclusion contour for tan/3 = 5, but only for 
A211 = 0.09, because the smaller couplings do not provide a 
region of 95% C.L. exclusion. 



FIG. 2. Exclusion contours at the 95% C.L. in the (mp,/mc, 
m x +) plane. The top figure is for tan/3 = 2, and three val- 
ues of A211, while the bottom figure is for tan/3 = 5 and 
A211 = 0.09. We give all contour plots as a function of the 
sneutrino mass. Because, for any given set of parameters, the 
sneutrino mass is very close to the smuon mass, the smuon 
contour plots lie very close to the sneutrino results, and are 
therefore not shown. 
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